Identifying Stellar Streams in the 1** RAVE Public Data Release 



R. KlementS B. Fuchs^ and H. W. Rix^ 
ABSTRACT 

We searched for and detected stellar streams or moving groups in the so- 
lar neighbourhood, using the data provided by the 1*** RAVE public data re- 
lease. This analysis is based on distances to RAVE stars estimated from a color- 
magnitude relation that was calibrated on Hipparcos stars. Our final sample 
consists of 7015 stars selected to be within 500 pc of the Sun and to have dis- 
tance errors better than 25%. Together with radial velocities from RAVE and 
proper motions from various data bases, there are estimates for all 6 phase-space 
coordinates of the stars in the sample. We characterize the orbits of these stars 
through suitable proxies for their angular momentum and eccentricity, and com- 
pare the observed distribution to the expectations from a smooth distribution. 
On this basis we identify at least four "phase space overdensities" of stars on 
very similar orbits in the Solar neighbourhood. We estimate the statistical sig- 
nificance of these overdensities by Monte Carlo simulations. Three of them have 
been identified previously: the Sirius and Hercules moving group and a stream 
found independently in 2006 by Arifyanto and Fuchs and Helmi et al. In addi- 
tion, we have found a new stream candidate on a quite radial orbit, suggesting an 
origin external to the Milky Way's disk. Also, there is evidence for the Arcturus 
stream and the Hyades- Pleiades moving group in the sample. The detections are 
further supported by analysing the stellar distribution in velocity and angular 
momentum space using the same Monte Carlo simulations. We also find that the 
significance of overdensities is comparable, independent of the space in which the 
stream search is conducted. This analysis, using only a minute fraction of the fi- 
nal RAVE data set, shows the power of this experiment to probe the phase-space 
substructure of stars around the Sun. 

Subject headings: Galaxy: solar neighbourhood — Galaxy: kinematics and dy- 
namics 
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Introduction 



According to current ideas about the cosmogony of galaxies these were assembled 
through hierarchical merging, which should result in a richly structured phase-space dis- 
tribution of dark matter and stars. Direct emperical evidence for such events has been 
sought among the stellar populations of the Milky Way with large scale surveys. Ongo- 
ing satellite accretion events hav e been discovered in several instances. A prominent ex- 
ample is the Sagittarius galaxy (llbata et al.l [l99J) with its associated tidal stream whic h 
wraps around the Ga l axy nearly perpendicula r to the galactic plane (llbata et al.l (120011 ) . 
Majewski et al.l (l2003r). iBelokurov et al. I (120061)'). Simm ilarly the Monoceros stream at low 



galactic latitudes (lYannv et ah 



2003 



Ibata et al.l 120031 ) or the recently discovered Orphan 



stream ( Belokurov et al. Hbood ) are i nterpreted as tidal debris from the Canis Major and Ursa 
Major II dwarf galaxies, respectively ( Penarrubia et alTlboOS : Fellhauer et al.|[2007 ). Numer- 
ical simulations have shown that such debris streams can survive as coherent structures over 



gigayears (jHelmi et al. I l2003l : iPefiarrubia et al. I l2005l : iLaw. Johnston &: Majewski 1 12005| ) 



Star streams, i.e. goups of stars on essentially the same orbits in the galactic potential, 
have also been detected as overdensities in the phase space distribution of stars in t he Solar 



neigh bourhood. The concept of "moving groups" originates from the work of Eggen (lEggen 



19961 . and references therein). Some of the moving groups are associated with young open 
clusters and can be naturally explained as clouds of former members, now unbound and 
drifting away from the clusters; these moving groups only reflect the nature of star formation, 
not necessarily that of hierarchical galaxy formation. 

However, data from the last decade seem to support the concept of co l d star streams in 



the Solar neighbourhood, consisting of old stars (5-10 Gyrs). iHelmi et al.l (119991 ) discovered 
the signature of a cold stream in the velocity distribution o f halo stars whi c h had been 
constructed fro m Hipparcos data. This was later confirmed by iGhiba fc Beers I (120001 ) using 
their own data. iHelmi et al.l (119991 ) argued that this stream formed part of the tidal debris of 
a disrupted sate l lite ga laxy accreted by the Milky Way, which ended up in the halo. Moreover, 
Navarro et al. I (120041 ) interpreted Eggen's Arcturus stream as such another debris stream, 
but in the thick disk of the Milky Way, dating back to an accretion event 5 to 8 Gyrs ago. 

Moving groups of old stars are also observed in the velocity distribution of thin disk stars 
in the Solar vicinity. Based on Hipparcos parallaxes and proper motions, iDehnen I (119981 ) 



found with statistical methods new evidence for the Sirius-UMa, Pleiades-Hyades, and Her- 
cules star streams. Even more convincingly these streams show up in the three-dimensional 



kinematical data of F and G stars in the Copenhagen- Geneva survey (INordstrom et al. 



20041 ). These moving groups are very probably not related to tidal debris streams, but orig- 
inated from dynamical effects within the disk itself like resonances with the inner bar of the 
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Milky Way fiDehnen II2OOOI ) or with spiral denisty waves (iQuillen fc Minchev II2005I ). 



Yet another str eam has been disco vered independently by lArifyanto &: Fuchs I (120061 . 
hereafter AF06) and iHelmi et al.l (120061 ) analyzing different data samples and applying dif- 
ferent detection techniques. The metallicities of the stars in this stream cover a broad range, 
from metallicities typical for old thin disk stars to that of thick disk stars. Thus the exact 
nature of this stream seems to be at present not fully understood. 



With the first data release from the RAVE collaboration (jSteinmetz et al.l 120061 ) a new 
large data sample of stars in the Milky Way became available, which is ideally suited for 
kinematical studies. To explore RAVE's potential, we analyze the velocity distribution of 
more than 7000 stars within a dis tance of 500 pc and search for overdensities in phase space 
using the projection technique of lAF06l . Even the first data release sample is so substantive 
that we can determine signal-to-noise ratios for orbital 'overdensities' and show the statistical 
significance of the detected overdensities. 



Data 



The RAVE DRl (Isteinmetz et all UQm) contains 25,274 radial velocities for 24,748 
individual stars together with proper motions and photometry from other major catalogs 
(Starnet 2.0, Tycho-2, SuperCOSMOS, USNO-B, DENIS and 2MASS). Its total sky coverage 
is ~ 4760 deg^. The only missing parameter for getting all velocity and position components 
is a distance estimate. The photometric data, however, allow the application of a photometric 
parallax relatioiil]. To calibrate such a relation for RAV E stars, we u sed main sequence (MS) 
stars (luminosity class V) from the Hipparcos catalog ( ESA 19971 ) with accurately known 
parallaxes (^ < 0.1) that could also be identified in the Tycho-2 USNO-B and 2MASS 
data. These other catalogs contain Bl, Vy and H-band magnitudecl. 



The absolute magnitudes of the stars were calculated through 

M = m - 10 + 5 log ■ ^ 



mas 



^We assume that the vast majority of RAVE stars are main sequence stars, because giants in the RAVE 
m agnitude ran g e wou ld he so far away, that the star density is substantially decreased. See also the discussion 



Siegel et all (120021 ) 



^We do not consider very metal-poor (sub)dwarfs, which can lie up to > 1 magnitude below the Solar- 
metallicity MS; for metal-poor stars we would overestimate the distance, and hence velocity. The effect for 
disk stars, however, will be negligible, because they exhibit approximately Solar metallicity. We will come 
back to the subdwarf problem later in the context of the detection of a new halo stream in the RAVE data. 
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We then considered two color-magnitude relations: [Vr — H) vs. My^ and [Bl — H) vs. 
Mbi- Figure [T] shows the relation Mvj, vs. {Vt — H). The dash-dotted line was chosen to 
remove all stars which seemed to be mis-classified as MS stars in the Hipparcos catalog. This 
cut was done by eye, since one could only reliably distinguish a (sub-)giant from a MS star in 
the vicinity of the MS by measuring its surface gravity. Next, the color-magnitude diagram 
was divided into three color bins in which we seperately fitted a color-magnitude relation 
(solid lines). For each bin, the intrinsic scatter of the color- magnitude relation was calculated 
as the 1/2 of the central 68% of the cumulative distribution of the differences AM between 
the true absolute magnitudes (obtained through the parallax) and those obtained through 
the fit. The scatter in each bin is denoted in Figure [H For our analysis we adopted the My^, 
vs. {Vt — -f/)-relation, because its intrinsic scatter is slightly smaller than that of Mbi vs. 
{Bl — H). Also, the errors of Bl for each star in the RAVE catalog are unknown, while the 
errors in Vr are given. The formal error in H (from 2MASS) is typically very small, ~ 0.01 
mag. For that reasons we chose the Mvj,{Vt — -ff)-relation to calculate absolute magnitudes 
Mvj, for the RAVE stars. If Vr was not given for a RAVE star, it was possible for most 
of the stars to get this magnitude through conversion formulas from other magnitudes like 
the USNO-B Bl and -Rl (Klement, unpublished). Figure [2] shows that a large fraction of 
all DRl stars has absolute magnitudes fainter than My^ = 6, i.e. lie in the section of the 
color-magnitude relation where its intrinsic scatter is smallest. Therefore we expect to get 
a large sample of stars with distance estimates good to 10% - 20% in most cases (see below 
and Figure [21, where the Myj,-distribution of our final sample is shown). 

We established a right-handed Cartesian coordinate system {x, y, z) centered on the 
Local Standard of Rest (LSR) with the x-axis pointing in the direction of the Galactic center, 
the y-axis pointing in the direction of Galactic rotation and the z-axis in the direction of the 
North Galactic Pole, (t/, V, W) denote the velocity components of a star in this Cartesian 
coordinate system. The Sun is assumed to lie at a distance Rq = 8kpc from the Galactic 
center and to move with velo city {U q,Vq,Wq) = {10. 0,5. 2, 7. 2)km with respect to the 



LSR (jPehnen fc Binney Ill998l ). A star's distance d follows from the apparent and absolute 



magnitudes {Vt,Mvj,) via Equation [H with its error ad determined through: 



1 



ad = -d\nlO^{av,)^ + {aMy^y, (2) 

where avj. and (Tmv^ denote the errors of Vr and Mvj, respectively. The main error is 
introduced by the intrinsic scatter of the color-magnitude relations (Figure [1]), which is 
translated into ctmvj.- 

The velocity components of a star can now be calculated from its position, radial 
velocity Vrad and proper motions {fia, fi's) through formulae given e.g. in Section III of 
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Johnson &: SoderblomI (119871 ): 
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(3) 



If /Iq, and are taken in units of [mas yr~^] and d in [kpc], the fac tor 4.74 gives the result 
i n [km s~^]. The coordinate transformation matrix B is defined in iJohnson fc Soderblom 
( 119871 ) and we evaluated it for the epoch 2000. The uncertainties in the ve locity components 
depen d on the uncertainties in radial velocity, proper motions and distance (jjohnson fc Soderblom 
1987h : 
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{4.7Mr[al + {f^^a,/dr] 
(4.74^)2 [< + (/.,a,/rf)2] 

bi2 ■ bi3 

^22 • ^23 
^32 ■ ^33 

The elements of the matrix C are the squares of the corresponding elements of the matrix 
B, i.e. Cij = bfj for all About 80% of the stars in DRl have radial velocity accuracies 

better than 3.4 km s~^ and 69.1% have a mean proper motion error of at most 2.6 mas yr~^ 
( jSteinmetz et al.l l2006l ) . If we restrict our sample to stars within dmax = 500pc from the 
Sun, i.e. within a volume of Vmax = ^d^ax, we can assume a typical distance for a star as 
the distance where the volume is half as large: < d >= O.b^^^dmax ~ 400 pc. If we further 
restrict the sample to stars with relative distance errors of ^ < 0.25 and assume a typical 
proper motion of 15 mas yr~^, we get as an estimate for the uncertainty of the transverse 
velocity <8.6 km s^^. The velocity error distribution in U and V is shown in Figure [31 It is 
peaked at an error smaller than 5 km s~^. 



For the kinematic stream search we selected only those stars which satisfied the following 
criteria: ^ < 0.25, d < 500pc, a total space velocity of \vtot\ < 350 km s~^, and au, cry < 35 
km s~^. We also restricted ourselves to stars with I > 200, b > 20, because in this field 
the observation density was highest. This leaves 7015 stars, among which we searched for 
signatures of local stellar streams. 



3. Search strategy for streams 



The distribution of stellar streams in veloc ity space is not trivial and depend s on the 
origin and the age, or type, of the stream stars (jPehnen Ill998l : iFamaey et al.ll2004j ). Princi- 
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pally, one has to distinguish between two types of streams: dynamical streams, i.e. groups 
of stars that are trapped in a small region of phase space by dynamical resonances, and 
streams, which stars originated from the same bound object, like a cluster or a satellite 
galaxy. Once those stars become unbound, they will have slightly different orbits, and hence 
frequencies, so that they will phase-mix. This leads to a broadening in the velocity dis- 
tribution, which is more prominent in the M^-component, because the vertical frequency is 
shorter than the horizontal frequencies. These streams will show a t ypical "banana" s haped 
distribution in U and V and a symmetric distribution in U and W (IHelmi et al.ll2006l ). The 
distribution in V of nearby stream stars with similar angular momenta is narrow, because 
these stars have basically the same azimuth. On the other hand, dynamical streams show 
a more clumplike structure in U and V. Their location is primarily set by the tangential 
velocity component V, because \^ is a measure for the guiding radius of a star in the Solar 
neighbourhood, whi ch itself is a measure for the location of Lindblad reso nances with a pe- 
riodic perturbation (IBinney fc Tremaine 1119871 : lOuillen fc Minchev II2005I ). They also show 
a broad range of W velocities due to phase- mixing (jPehnen I Il998l ). Generally, scattering 
processes and measurement errors tend to symmetrize and broaden the velocity distribution 
of stellar streams of both types and makes it difficult to identify them in velocity space, 
particularly if the stream is very widespread spatially and consists only of a small fraction 
of stars compared to the smooth background. 



According to iHelmi et al.l (119991 ) and iHelmi fc de Zeeuw I (120001 ) a better approach to 
search for streams is in the space spanned up by the integrals of motion - i.e. the en- 
ergy E and angular momenta L^, L± = ^/L^ + or Ltotai- For example, iHelmi et al. 



(1l999l ) detected two fossil streams originating from the same progenitor as a clump of so- 
lar neighbourhood halo s tars i n {Lz,L±) space. Using cosmological N-body simulations, 
Choi. Weinberg fc Katz I (120071 ) showed that over 8 Gyr the overall position of a disrupting 
satellite and its tidal tails basically remains the same in (E, Lz) space, but the shape of 
the distribution shifts and one satellite can produce several apparently disassociated clumps. 
This, together with the fact that the total angular momentum is not really an integral of 
motion, can obscur e the signature of a well defined moving group in phase space. Therefore, 
Helmi et al.l (120061 ) proposed to look for stellar streams in a space spanned up by apocen- 
ter, pericenter and angular momentum Lz. Moving groups cluster around lines of constant 
eccentricity. 

We follow a similar approach in the sense that we try to find stars with the same orbital 
eccentricity. Our strategy for finding nearby stellar st reams in veloc ity space is based on 
the K eplerian approximation for orbits developed by iDekker I (119761 ) and summarized in 
AF06I . We assume an axisymmetric potential and that stars in the same stellar stream move 
on orbits that stay close together, which is justified by numerical simulations of satellite 
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disruption (iHelmi et al.ll2006l ). These stars should form a clump in the projection of velocity 
space spanned up by ^/lP~+'2V^ and V. The latter is related to the angular momentum 
which defines the guiding center orbits of the stars. The first quantity is a measure for a 
star's eccentricity e. For a flat rotation curve e is given by 



— • 

Vlsr denotes the circular velocity of the Local Standard of Rest for which we take the lAU 
standard value 220 km s~^. 

Also, it can be shown that the radial action integral Jr in the Keplerian approximation 
is approximately equal to ^^^{U"^ + 2V^^) (Fuchs, private communication). Therefore, the 
quantity \^LP + 2V^ should be robust against slow changes in the potential. 

For the current investigation we do not include W velocity components in our search 
for local streams, because - as mentioned above - strong phase-mixing will smear out any 
coherent feature over short timescales. However, in the appendix Section [5] we will make 
use of the W velocities when we compare our method to the more traditional searches in 
{U,V,W) and {Lz,L±) space. 



4. Results and Discussion 

We now use the distribution of the 7015 selected RAVE stars in \^TP~+'2V^ vs. V, 
shown in Figure HI to search for substructure in the kinematic distribution. Figure H] also 
shows the error ellipses for some points; most of the relative errors are rather small and even 
though some of the larger errors produce an uncertainty in the exact position of the data 
points and tend to smear out substructure, several suggestive " overdensities" of stars are 
visible by eye. In this section we focus on identifying overdensities and on quantifiying their 
significance. 



4.1. The Wavelet Transform 



To this end, we follow the same procedure outlined in IAFOGI and use the wavelet- 
transform technique using a two-dimensional analysing wavelet ^'(x, ?/). We bin the data 
in pixels of 2 km s~^ width on each side and calculate the value of the wavelet-transform in 
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each bin through: 

w{x,y) ~ J J dx'dy''^{x — x' ,y — y') X 



N-l 

^5{x' - Xi)5{y' - yi 



i=0 
N-l 



^■^{x - Xi,y - yi), (6) 



i=0 

where = 7015 is the number of stars in our sample. 



Motivated by the work of lSkuljan et al.l (119991 ). lAF06l used a Mexican hat shaped kernel 



function to detect overdensities in y/W^ + 2V^ vs. V. However, since the former quantity 
is not uncorrelated to the latter, we expect clumps that would be roughly spherical in, e.g., 
U vs. V, to be elongated along U = 0. That means that the performance of a Mexican 
hat shaped kernel is sub-optimal and might influence the significance of the overdensities. 
Therefore, we build a kernel function that is elongated about a factor q along the V-shaped 
U = lines in the following way: 

We rotate the coordinate axes such that the new \^-axis, V, lies along the line U = 
and the new \^U"^ + 2V'^ axis perpendicular to V. For simplicity we rename V as x and 
\UJ^+2V^ as y. This implies two different rotations, depending on wether the features 
are elongated in the region ^ < (clockwise rotation by the angle (f = arccos or in 
the region > (counter-clockwise rotation by (f). Based on the Mexican hat shaped 
function, we then express the analysing wavelet as a function of the new coordinates {x',y') 
and elongate it along the x'-axis by a factor q: 

*<^'f''-(2-w-5)^^K-2(^-|^)- p) 

This function is normalized in the sense that it's volume integral is zero. The scale parameter 
a is a measure for the extend of the "bumps". With the transformation equations 

x' = \ T \l^y (8a) 
y' = ±y|x+y|y, (8b) 
where the upper sign stands for the case V < 0, it's straightforward to calculate the kernel 
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in the unrotated coordinate system: 

^{x,y) = (2-^^{{l + 2q^)x^ + {2 + q^)y'T{l-q')2V2xy) 
exp(-^^((l + 2q^)x' + (2 + q^)y^ ^ (1 - q^)2V2xy) 



X 



(9) 



Then Equation ([6]) can be used to calculate the value of the wavelet transform in each bin. 

For the elongation parameter q we used -\/3, because the projection of a certain range 
in V along the line t/ = is \/3 times as long. For the scale parameter a we chose 7 km s~^, 
comparable to the mean errors, which extracts the overdensities most clearly (Figure [5])|f|. 

Figure [5] shows that the bulk of our sample stars have thin-disk-like kinematics, with \ U\ 
and \ V\ < 20 km s~^ The contour levels do not reflect some finer structures at —80 ^V<20 
km s~^. In this range we would expect the Hercules (V ~ —50 km s~^) as well as the Hyades- 
Pleiades {V ~ —20 km s~^) and Ursa Major/Sirius {V ~ +4 km s^^) streams (IFamaey et al. 



2004J ). Indeed, the distribution in {y/lP + 2V^, V) is somewhat bulged at these velocities. 



Among the stars with high orbital eccentricities, presumably thick disk or halo stars. 
Figure [5] appears to show five different overdensities. The first on e at V ~ —95 km s~^ 



i s mos t probably the stream which was discovered independently by IAFOGI and iHelmi et al. 
(120061 ). The second overdensity at V — 120 km s~^ has the same kinematics as the Arcturus 



group discovered by iNavarro et al. I (120041 ) . These two streams have most probably an origin 
external to the Milky Way disk, be cause their phase- space distribution resembles that of 
simulated accreted dwarf sattelites (IHelmi et al.ll2006l ). However, nothing is known about 
their progenitors yet. The third feature at ^ ~ —160 km s~^ has to our knowledge not yet 
been described in the literature. The same is true for the two clumps at y ~ —200 km s~^ 
and V ~ +50 km s~^. However, at these velocities our sample is very sparsely populated, 
so that just a few stars are enough to create a high value of the wavelet transform. Every 
isolated star in the middle of a bin increases the value of the wavelet in this bin by 2. This 
also explains the other clumps lying off the V-shaped main feature in Figure [5l 



4.2. Subtracting a Smooth Velocity Distribution 

To test whether any of these kinematic overdensities, or streams, reflected as peaks 
or clumps in the wavelet transform, are significant, we performed 250 Monte Carlo (MC) 



■^We also tested different choices for the set {q,a), hkc (\/2, 8km s ^) and basicaUy got the same results. 
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simulations of the same number of stars as in our sample, wh ich we randomly draw from a 
Galactic model consisting of three Schwarzschild distributions (IBinney fc Tremaine 1119871 ) to 
represent the thin and thick disk and the halo. The goal was to create a "smooth" reference 
model velocity distribution that matches the overall velocities of the RAVE sample. To this 
distribution, we then added normally disributed velocity errors, based on the observed error 
distributions shown in Figure [31 We chose a local thick-to-thin disk and halo -to-thin disk 



norma li zation of 0. 1 and 0.001, respectively, in agreeme nt with the v a lue fr om lChen et al. 



( I2OOII ): lJuric et al.l (120081 ) and the upper limit given by ISiegel et al.l (120021 ). The thin and 
thick disk and halo are assumed to have velocity dispersions {au, cry, (Tw) and rotational 
offsets from the LSR equal to (25,21,17,-5), (74,50,50,-44) and (189,97,100,-219) km s-\ 
respectively. The values for the thin disk have been chosen to best match the smooth part 
of the observed velocity distri bution (Figure El), althoug h the ratio oujoy in the thin disk 
should be approximately 1.6 (jPehnen fc Binney I Il998l ). Also, we chose an offset in the 
VT-distribution of +3 and +10 km s~^ for the thin and thick disk to get a good match0. 
The matches shown in Figure [61 are not perfect, but the slope of the velocity distribution 
at negative V velocities is reproduced well in the simulation. The distributions differ most 
clearly around V = km s~^, because - possibly due to the Sirius moving group - the RAVE 
stars are peaked around V ~ +5 km s^^. Moreover, the assumption of a Schwarzschild 
distribution does not reflect the skewness of the velocity distribution due to the asymmetric 
drift effect. 

For each MC sample of 7015 stars we binned the velocities and calculated the wavelet 
transform in each bin the same way as for the real sample. Figure [3 displays one randomly 
chosen MC realization. One can clearly detect several peaks in the wavelet transform, seem- 
ing "overdensities", which must be, however, by construction due to Poisson noise. This 
shows that also in our real sample the possibility for such "fake" streams exists. 

We proceeded with calculating the mean value Wi^j and the standard deviation cTjj- of 
the 250 MC wavelet transforms in each bin (i, j). Because many bins are not populated with 
stars, the standard deviation in those formally has values cxjj < 1; when this was the case, 
we set (Tjj- = 1. The contours of the mean and the standard deviation are shown in Figure [HI 
and Figure [9], respectively. 

While a single Monte Carlo sample shows v arious kinds of clumps and overdensities due 
to statistical fluctuations (Figure [TJ and see also lAF06l ) . the mean value of all wavelet trans- 
forms represents a very smooth distribution, because the clumps disappear when averaging 



''This deviation from the expected mean value of 14^ = can be due to the hmited sky coverage of our 
sample 
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over all samples. 

The standard deviations aij define the significance of structure in our real sample in 
the following way: Because we are only interested in overdense regions, i.e. pixels where the 
wavelet transform takes on positive values, we set the value of the wavelet transform of the 
data, w°^j^, as well as the mean value of the wavelet transform of the 500 MC samples, wfj'", 
to zero, whenever it is < 0. By doing so, we make sure, that the residuum — wf^j'") is 
> in each pixel, where our "smooth" model contains no stars. By dividing the residual 
(wg^ - wg^) in each pixel through cTjj- we define the significance of peaks of wg*. For a 
positive stream detection, we require a significance of at least 2. Figure [TO] shows the result. 



4.3. Significance of the Streams 



The first thing to mention is that in all bins where the standard deviation has a value 

. In these bins, one 



of 1, the significance has the same value as the residual (wg'^ 
or a few stars will result in a significant overdensity with cxjj > 2, since for every star in 
the middle of a bin, the value of increases by 2. This effect could then produce "fake" 
stellar streams, just because of Poisson noise in our RAVE sample and the smoothness of 
the combined Monte Carlo samples. Since a solution to this problem is difficult, we avoid 
speculating about the overdensities in bins with aij = 1, which we encircled red in Figure [TUl 
We only give the number of stars that are contained in this features, as derived from the 
scatterplot. Figure |H For example the extended clump at ~ —200 km s~^ corresponds 
to a group of 8 stars, while the one at ~ —230 km s~^ contains only 2 stars. Further 
investigations with more data could provide more hints for the nature of this features. 

We try to match the remaining statistically significant peaks with streams that have 
been already described in the literature. 

Centered at ^ ~ +4 km s~^, \U\ ^10 km s~^, there is a bi g clump, signif i cant at a level 



W, 



of> 10. This is probably the Ursa Major/Sirius moving group (jPehnen lll998l : lFamaey et al. 



20041 ). its signal possibly amplified through the inability of our simple three-component MC 
models to fully match the observed thin disk distribution (see also Figure [6]). The core of 
this group is in the direction of Ursa Major. It's members are distributed all around the 
sky and can be very close (Si rius at 2.65 pc distance is a member of the group), putting 
the sun inside of the group (IBannister fc Jameson 1 120071 . and references therein) . The 
current understanding is that the Sirius stream consists not only of a cluster of coeval stars, 
but also of different kinds of field stars, whi ch altogether could have been forced on similar 



orbits by a spiral wave gravitational field (jSellwood fc Binney 1 12002| ; iFamaey et al.l 12004 
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Quillen fc Minchev 1120051 ) . 



The elongated clump streching from V ~ —50 to ~ —75 km s^^ with a peak signif- 
icance of 0" = 4.5 can probably best be described as the Hercules stream. This stream can 
be explaine d as the scatte ri ng of stars off the G alactic bar induced by the Outer Lindblad 
Resonance (jPehnen Il2000l ). iFamaey et al.l (120041 ) also find a group of most likely thick disk 
and halo stars located at = —53.3 ± 41.36 km s~^. Further, at = —60 km s~^ there 
exists the moving gro up HR1614 which is thought to be a dispersed open cluster due to its 
chemical homogenity (lEggen Ill996l : lde Silva et al. 1120071 ). It is possible that these groups are 
present in our data, too, amplifying and elongating the signal of the Hercules stream. Better 
velocity estimates would be needed in order to clearly distinguish between these features. 



The feature at ^ —100 km s ^ stands out at the 4.3cr- 



stream discovered independently from lAF06l and iHelmi et al.l (120061 ). Most probably, this 



evel. It corresponds to the 



stream has its origin outside of the Galaxy and is a relic from the build-up of the Milky Way. 

At a velocity of ~ —160 km s~^ there is an overdensity which stands out at the 
3.0cr-level in the center. A comparison with Figure reveals that this feature is probably 
more elongated, like one would also expect for a moving group of stars at such velocities. 
However, with the current sample size, our method seems not to be fully able to recover 
the whole range from V = —180 km s~^ to ^ = —140 km s~^ as a statistically significant 
overdensity. Figure [TT] shows the stars that make up the new feature in a color-magnitude 
diagram consisting of — if color and absolute magnitude My^, . Stars that lie in the area 
with a significance greater than 2 are plotted with thick asterisks, while those that lie in the 
range V = —180 km s~^ to = —140 km s~^ and presumably also belong to the stream have 
thin asterisks. We also show isochrones for a 13 Gyr old population of stars with metallicities 
of [Fe/H]=-1.5, -1.0, -0.5 and 0.0 (from left to right). Most of the stars seem to be consistent 
with having Solar metallicity, however this is a result which would be expected from the 
way we calculated their distances and velocities. Our assumption was that all stars are MS 
stars of luminosity class V, and therefore we a priori excluded the possibility to have very 
metal-poor (sub)dwarfs with fainter luminosities (see also footnote 2 in Section [2]). Right 
now, we can only wait for future data releases providing metallicity estimates to solve this 
issue. 

In appendixRIwe in vestigate the Geneva-Copenhagen survey of the Solar neighbourhood 
(INordstrom et al. 1 12004| ) with our method; although we do not determine the significance 
of overdensi t ies w ith the help of MC sampling, we can see the same strong features that 
Helmi et al.l (120061 ) detected in their paper. It is striking that the same stellar streams can 
be detected in different datasets using different methods. 
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Therefore - through the detection of 3 aheady known stellar streams in our small RAVE 
dataset - we have enough confidence to claim that the new feature centered on ^ = — 160 km 
s~^ is a strong candidate for a new stellar stream. The IV-velocities of stars in this feature 
show a wide range of values, which makes it likely, that this stream was accreted long ago as 
tidal debris during the formation of the Galaxy. To draw more conclusions, one would need 
metallicity measurements, which will be available from RAVE in the future. 

It is possible to give constraints on the density contrast in the streams in the Milky 
Way stellar halo. Therefore we estimate how many stars are needed to give a significant 
stream detection and divide this number by the total number of halo stars. Because from 
kinematics alone it is very difficult to del imit halo from t hick d isk stars (and even with 



metallicities there is no clear boundary, see IChiba &: Beers I (120001 )). we simply take all stars 



outside of (V ± 2cry)^j^.^ ^.^j^ = —5 ± 42 as halo stars and accept a contamination from the 
thick disk. This gives 554 stars, of which 230 (42%) are part of the Hercules stream, 47 (9%) 
belong to the AF06 stream and 7 (1%) make up the new detected stream. We conclude that 
a few percent density contrast are enough to give a significant stream detection. 



5. Searching for Stellar Streams in {U,V,W) and {Lz,L±) space 



The t raditional con cept of moving groups refers to distinct clumps in velocity space 
{U,V,W) (lEggen Ill996l . and references therein). However, as explained in Section [3l this 
space is not immune to the effects of phase-mixing: over time, distinct clumps will dis- 
perse due to small initial differences in the oscillation periods of the constituting stars and 
encounters with e.g. molecular clouds. For a part of a stellar stream to be in the Solar 
neighbourhood at the same time means its stars must have similar V motions, but the U 
and W velocities can differ by hundreds of km s~^. 

In contrast, stellar streams are able to remain as cohere nt features in the space of 
integrals of motion. Therefore, in the style of Helmi's work (IHelmi et al.lll999l ). we will 
investigate to what extent the detected streams from Figure [TU] can be recovered in (f/, V, W) 
and {Lz, L±) space. 



5.1. Analysing the RAVE data in {U,V,W) space 

In Figure [12] we show the distribution of our 7015 RAVE stars in velocity space. The 
upper panels show how the stars are distributed with respect to their (f/, V, W) velocities. 
The distribution is far from smooth, as can be seen more clearly in the lower panels, where 
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we display the wavelet transform of the velocity components. For the convol ution we chose 



a Me xican hat kernel with a scale parameter of a = 5 km s ^ (similar to ISkuljan et al. 



19991 ). The contour levels range from 1 (black) to 50 (red). The most striking features are 
located aound V = 0, U <= and most probably can be attributed to dynamical thin disk 
streams like the Sirius moving group. The elongated shape along the [/-axis is typical for 
stellar steams and is caused by the fact that we observe stars with slightly different orbital 
phases. The classical moving groups Sirius at {U,V) ~ (+6, +4) km s~^, Hyades-Pleiades 
at {U, V) ^ (—25, —15) km s~^ and Hercules at (U, V) ^ (—2 , —50 ) km s"^ can be found 



with good agreement with the values given by iFamaey et al.l (120041 ). The incomplete sky 
coverage of our sample as well as the already described effects of phase mixing make it hard 
to speculate about the meaning of the many other clumps. We nevertheless performed a 
statistical Monte Carlo analysis analog to that described in Section 14.21 using 200 Monte 
Carlo simulations. Figure [T^ provides the significance map of the velocity space, showing all 
overdensities from Figure [T2l that have a significance of at least 2. 



Compared to (V, VU"^ + ^V"^) space (c.f. Figure [TO]) there are more features present, 
but most of them are very small and at lower significance levels. Besides the prominent 
red feature around {U,V) ~ (0,0), which can be attributed to the Sirius stream, it is not 
easy to associate the identified streams from Figure [TOl with clumps in the ([/, V) or {U, W) 
diagrams, respectively. Because old moving groups are not expected to remain coherent as 
small clumps in velocity space, we likely over-resolve substructure due to the small deviations 
of our Monte Carlo model from the data. These deviations can also lead to a wash-out of 



the b anana-shaped structures in [U, V) that are expected for old tidal streams (IHelmi et al. 



20061 ). Large errors in the velocities can wash out these structures, too, but in our case the 
errors are small enough and should allow their detection. The fact that a clear classification of 
particularly old stellar streams from the {U, V) and {U, W) distributions is not possible, and 
that many small features are even less significant than the overdensities in (V, y/U"^ + 2V'^) 
space lead us to conclude that (f/, V, W) space is not suited well to detect moving groups. 



5.2. Analysing the RAVE data in {Lz,L±) space 

In a spherical static potential the two components of the angular momentum and 
L± = (L^ + Ly)^ are integrals of motion. In this case, tidal and dynamical streams remain 
as coherent clumps in {Lz,L±) space. The asumption of a spherical potential seems valid 
for large heights above or below the disk which are reached only by halo stars, but not 
necessarily in our case of nearby stars with planar orbits. These stars move under the 
influence of a flattened potential and the radial action integral Jr would be better suited 
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as a conserved quantity. The space of integrals of m otion has been shown to be useful for 



the detection of stellar halo streams in recent studies (IHelmi et al.lll999l : iHelmi &: de Zeeuw 



2OOOI : IChiba fc Beers Il2000[ ). 



In Figure UM we show the distribution of our RAVE stars in {Lz,L±) space, again as 
a scatter plot (upper panel) and after wavelet transformation using a Mexican hat kernel 
with scale parameter a = 45 kpc km s~^. For stars near the sun — RoiV + Vlsr), so we 
expect to identify again the overdensities from Figure O 

A comparison shows first of all that [L^, L±) seems to be much more structured than 
{V, ^/U"^ + 2V'^) space, but this is partly due to the applied chosen scale parameter a in the 
kernel function. The other reason is that now all three velocity components go into the 
calculation of the angular momentum components. Also, there are no "forbidden regions" 
as in the case of (V, ^/lP~+'2V^) space. 



Some overdense regions in the lower panel can be related to overdensities in (V, ^/[P + 2V^) 
space. For example, the "bump" in the distribution in Figure |5] at V ~ —15 km s~^ is also 
present in Figure [TH at Lz ~ 1640 kpc km s~^. Also, the Hercules stream is present at 
Lz ~ 1350 kpc km s~^. 

We proceeded to investigate which overdensities proof significant in the same way as 
before by using 200 MC samples that were drawn from the 3-component Galaxy model 
described in Section 14.21 After building the mean value and standard deviation of their 
wavelet transforms we derived the significance map, which is shown in Figure [151 

As in the case of {U, V, W) space the number of features with significance a > 2 is larger 
than in Figure [10 We think this is caused by the extra-information of W velocity required 
for the MC samples. Some of these features are very small and again belong to regions of 
phase space that are sparsely sampled, but others seem to be comparable to the already 
known ones. We encircled all features which correspond to the detections of Figure [10] with 
black circles for the appointed streams or with red circles for the doubtful ones. Additionally, 
there is a new feature present at Lz ~ 1650 kpc km s^^ which we identify as the Hyades- 
Pleiades group. A hint on the presence of this group has already been seen as a "bulge" 
in the distribution of stars in {V, \^LP + 2V^) space (Figure [5]). It is interesting that some 
distinct features in {Lz,L±) space (e.g. the new stream) project onto the same feature in 
(V^, y/lP + 2V^) space. Furthermore, when we select stars from the fastest rotating clumps 
(Hyades-Pleiades, Sirius and the red encircled clump on the right side in Figure [TBI) we 
find that they span a wide range of \/U^ + 2V'^ values. On the other hand we also find 
that clumps in (V, y/lP + 2V^) space are not necessarily projected onto clumps in {Lz,L±) 
space, indicating that L± is not conserved for their motion. Therefore, we conclude that L± 
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compared to \/U'^ -\- 2V^ is not a good choice as an approximate integral of motion if the 
disk's potential dominates the movement of the stars. 

We note as a main result that the significance of the detected features is comparable, 
no matter which space we use for the stream search. As explained before, we expect the 
significance of the detections to go up with a larger sample size. For best results, it is probably 

good to combine the searches in {Lz,L±) and (V, \/lP~+'2V^) space to find stellar streams 
on both inclined and planar orbits. Also, the significance is influenced by the statistical 
method which is used to identify overdense regions. Our MC method has the drawback that 
a good Galactic velocity distribution model is needed in order to avoid " fake" detections in 
regions where the data points are rare. Maybe a combination of different analyses can also 
enhance the detection efficiency. 

6. Conclusions 

To explore RAVE's practical potential for finding kinematic substructures (streams, or 
moving groups) in the Solar neighbourhood, we have studied a sample of 7015 nearby stars 
from the RAVE first public data release, which was selected to have acceptable distance 
and velocity estimates. Distances have been derived from a photometric parallax relation, 
based on Hipparcos MS stars using the Vr and H bands. For our sample, we required a 
distance error better than 25%, while the mean velocity errors in U and V are ~ —5 km s~^, 
respectively. 

To search for streams, we plotted V velocities against the quantity y/W^ + 2V^, which 
essentially corresponds to angular momentum Lz against orbital eccentricity e. Streams 
should form clumps in this projection of phase space. To identify such clumps more obvious, 
we applied a wavelet transform with a modified two-dimensional Mexican kernel as the 
analysing wavelet. Several overdensities are visible in this sample, presumably corresponding 
to stellar streams. We tested whether these clumps are real or due to Poisson noise using 
250 MC simulations drawn from a three component Schwarzschild distribution. While it 
is very difficult to draw conclusions about the amount and significance of sub-structure 
present in the very few halo stars, we detected three already known stellar streams and one 
candidate for a newly discovered stream. The latter is present as a broad feature in the range 
—180 <V< —140, centered a,t V ^ —160 km s~^ and from its kinematics would belong to 
the stellar halo population. Its other velocity components (high W velocities with a mean 
of W = 117 ± 26 km s^^, U = 6.7 ± 55 km s~^) make it likely that this stream is part of the 
tidal debris from an accreted satellite rather than a dynamical resonance. The kinematics, 
however, could be somewhat biased, because we did not correct for metallicity effects in the 
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photometric parallax relation. The other moving groups are the Sirius s tream at V ^ +4 km 



s ^ , the Hercules stream cent ered a t V y —65 km s ^ (IDehnen I Il998l : iFamaey et al.l 120041 ) , 
and the stream discovered by AF06 and Helmi et al. ( 20061 ) atV ~ —100 km s~^. 



We are missing a detection of the Hyades-Pleiades and Arcturus streams which were 
still present as peaks in the wavelet transform of the data in Figure [51 Here, the difference 
— wf^j^) is too small to make the peak in wf^^^ statistically significant. 

However, we were able to detect the Hyades-Pleiades moving group in velocity and an- 
gular momentum space, which we used to complement our stream search. By comparing the 
detections in {U,V,W), {Lz,L±) and (y, vf/M-2T^) space we showed that i) the signifi- 
cances of the features were comparable, ii) that velocity space is sub-optimal for detecting 
stellar streams, but may be useful as an additional information on the age and origin of a 
stream (tidal streams look differently from dynamical streams, old streams disperse) and iii) 
that L± is a poor approximate integral of motion compared to yLP + 2V^ if the stars move 
under the main influence of the disk's potential. 

The fact that only a fraction of the DRl stars were enough to find significant sub- 
structure in the Solar neighbourhood shows the power of the method we used. Subsequent 
data releases will not only enlarge the sample of stars and eventually allow the detection 
of further clumps in phase space, but also provide measurements of log (7, Tpjf and \M/H]. 



The full dataset will include these measures for up to one million stars (jSteinmetz et al. 



20061 ). Simply speaking, and neglecting an increase in the spatial coverage, an enlargement 
of the sample size by a factor of N leads to N times more stars in a stellar stream, while the 
standard deviation of the MC samples with be increased by a factor of (Poisson noise). 
The mean value of the wavelet transform of the MC samples stays the same, so that the 
significance of the features will be N/ \/N = ^/N times as high. In this way we can hope 
to detect more features and put further constraints on the exact nature and origin of the 
detected streams. 



A. Streams in the Geneva-Copenhagen survey 



We show the performance of our method on the data take n from the Gen e va-Co penhagen 



survey of the Solar neighbourhood ( iNordstrom et al. I I2004J ) . iHelmi et al.l (l2006l ) searched 



for stellar steams in velocity and APL space in these data. While the dynamical streams 
(Sirius, Hyades-Pleiades and Hercules) can be seen as clumps in velocity space, a detailed 
statistical analysis of the APL space revealed arou nd 10 further overd ensities oriented along 
two to three segments of constant excentricity. iHelmi et al.l ( 120061 ) divided the stars in 
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the overdense regions into three groups based on their metalhcities. There seems to be 
a relation between these three groups, the Arcturus stream and the stream detected by 



Arifvanto fc Fuchs I fl2006f ). 



In Figure [16] we just show the distribution of the 13240 Geneva-Copenhagen stars in 
V vs. y/U"^ + 2V'^, after convolution with the analysing wavelet given by Equation ©El 



The most striking features are those that have also been found by lHelmi et al.l (l2006l ). Ad- 
ditionally, a comparison with our RAVE sample (Figure seems to suggest that there is 
much more sub-structure present in the Geneva-Copenhagen survey, although it contains 
only ~ 2 times more stars. We find that 10582 of the 13240 stars are included in these 
features with a value of the wavelet transform > 1. That makes up for 79.9% of the 
sample. Our RAVE sample contains 7015 stars, of which 5601, or 79.8%, constitute to re- 
gions with Wij > 1. So it seems that the size of the sample is the main reason that in the 
Geneva-Copenhagen survey we find more distinct overdensities. Furthermore, if we relate 
the number of patches above the V-shaped U = line (30 in the Geneva-Copenhagen, 14 
in the RAVE sample) to the number of stars in the sample, we find a similar relation; the 
larger number of stars increases the probability to find outliers in less populated regions. 
Finally, stars in the Geneva-Copenhagen survey have accurate trigonometric parallaxes from 
Hipparcos with relative errors (ct^/tt) better than 10 per cent, in contrast to the less acurate 
photometric parallaxes of our RAVE stars. This implies less accurate velocities that tend to 
smear out small nearby features in phase space. 

Although it's unlikely that the majority of the small features in Figure [16] has any sta- 
tistical significance, this example still shows that a large sample size together with good 
distance estimates is crucial for detecting stellar streams in phase space. Still, our photo- 
metric parallaxes seem to be good enough and our sample large enough for significant stream 
detections. Without investigating the statistical significance of the features in Figure [161 we 
leave this example as further evidence that stellar streams will be detectable as clumps in 
{V, + 2V^) space. 
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Fig. 1. — Absolute magnitude calibration for Hipparcos MS stars with a parallax accuracy 
better than 10% with Vr — as the color. All stars above the dash-dotted hne are assumed 
to be giants incorrectly classified as MS stars and are not considered for the color-magnitude 
fit. The solid line shows the My^ — (Vt — H) relation adopted, and the dispersions a about 
this mean relation are indicated in three color regimes. 
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Fig. 2. — Distribution of absolute magnitudes of all DRl RAVE stars (solid line) and our 
selected sample (dash-dotted line). A large fraction has Mvj. > 6. 
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Fig. 3. — Distribution of errors for the velocity components U (dashed) and V(dotted) for all 
stars in RAVE DRl as derived from Equation |H These errors include radial velocity errors, 
proper motion errors and distance uncertainties. Both distributions have a peak around 2-3 
km s~^ and a long tail reaching to much higher values. 



Fig. 4. — Distribution of our sample of RAVE stars in \/TP~+'2V^ vs. V. We also show 
error ellipses for a small subset of stars; note that while the large error ellipses are most 
prominent, most ellipses are very small. 
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Fig. 5. — Contours of the wavelet transform of the distribution of our sample of RAVE stars 
(Equation [6]) in y/lP + 2V^ vs. V. The contour levels are displayed in the color bar, the 
scale parameter of the analysing wavelet is 10 km s~^. 
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Fig. 6. — Velocity distributions of our selected RAVE stars compared to one MC realization. 
The MC sample consists of three Schwarzschild distributions for the thin and thick disk and 
the halo, respectively. Note the peak in the data at around V = +4 km s~^, which we later 
interpret as the Sirius moving group. 
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Fig. 7. — Same as Figure [5l but now for one of our 250 Monte Carlo samples of 7015 stars 
drawn from a smooth velocity model (Section |4]). Due to Poisson noise, some overdensities 
emerge. 
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Fig. 8. — Same as Figure Ui but now the mean value of the wavelet transform for all 250 
Monte Carlo samples is shown. Only values > 1 are displayed. 
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Fig. 9. — Same coordinates as Figure Hand [8], but here the standard deviation of the wavelet 
transform among the 250 Monte Carlo samples is shown. Outside of the contoured region, 
the variance is set to unity. This variance map is used to asses the significance of wavelet 
transform peaks in Figure [51 
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Fig. 10. — Significance of tlie overdensities seen in Figure [5], obtained as described in tfie 
text. Note tliat only areas witli a > 2 are displayed. We encircled black and labeled all 
features which we considere to be stellar streams. We are not making further suggestions 
about the four features encircled with red, but instead give the number of stars which make 
out these features. 



-SO- 




LO 1.5 2.0 2.5 3.0 3.5 4.0 

Vt-H [mag] 



Fig. 11. — CMD of the stars that make up the feature between V ^ —140 and V ^ —180, 

highUghting those in the significantly overdense region —167 km s~^ < V < —153 km s~^ 
(thick asterisks). Also, isochroncs for a 13 Gyr old population with mctallicity [Fc/H]=— 1.5, 
— 1.0, —0.5 and 0.0 are plotted. The linearity of our adopted photometric parallax relation 
in the range Vt — H > 1.9 is clearly visible. 



Fig. 12. — Distribution of RAVE stars in velocity space. The upper panels show scatter 
plots of U vs. V and U vs. W, respectively. A lot of sub-structure can be seen. In the lower 
panels we show contours of a wavelet transform of the data. The convolution has been made 
using a Mexican hat kernel with scale parameter a = 5 km s~^. 
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Fig. 13. — Significance of tlie overdensities seen in Figure [121 Only areas with cr > 2 are 
displayed. Note the richness of features compared to Figure dUl with many having a very 
small extension. 



Fig. 14. — Distribution of RAVE stars in {Lz,Lj_) space. The upper panel is a scatter plot 
of the data. In the lower panel we show contours of a wavelet transform of the data. The 
convolution has been made using a Mexican hat kernel with scale parameter a = 45 km s~^. 
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Fig. 15. — Significance of tlie overdensities seen in Figure [TH Only areas witli cr > 2 are 
displayed. Black circles mark features which we connect to the detections of Figure [TU| there 
is further evidence for the now significant Hyades-Pleiades group. Red circles mark features 
which can be connected to doubtful features in Figure [TOl Note that for stars near the sun 
Lz ~ RqV , were we take Rq = 8 kpc. 
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Fig. 16. — Contours of the wavelet transform of 13440 stars from the Geneva-Copenhagen 
survey plotted in V vs. VU"^ + 2V'^ (most of wh i ch ar e not statistically significant). The 
features that have been described by iHelmi et al.l (120061 ) have been labeled. 
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